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Abstract 

A quantitative description of the configurational part of the total energy 
of metallic alloys with substantial atomic size difference cannot be achieved 
in the atomic sphere approximation: It needs to be corrected at least for 
the multipole moment interactions in the Madelung part of the one-electron 
potential and energy. In the case of a random alloy such interactions can 
be accounted for only by lifting the atomic sphere and single-site approx- 
imations, in order to include the polarization due to local environment 
effects. Nevertheless a simple parameterization of the screened Coulomb 
interactions for the ordinary single-site methods, including the generalized 
perturbation method, is still possible. We obtained such a parameterization 
for bulk and surface NiPt alloys, which allows one to obtain quantitatively 
accurate effective interactions in this system. 
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I. INTRODUCTION 



One of the main problems of modern alloy theory is to establish a quantitatively accu- 
rate description of the configuration-dependent part of the free energy, i.e., the difference 
in the total energies of alloys with different atomic arrangements on the underlying lat- 
tice, in terms of effective cluster interactions which may subsequently be used in statistical 
thermodynamics simulations 0^0] • Even without lattice relaxation effects (which are not 
considered here, although they play an important role in the phase equilibria of many alloy 
systems) a solution to the problem is still a challenge especially in the case of inhomoge- 
neous systems such as surfaces in the presence of long-range and multi-site interactions 
which cannot be neglected. 

The challenge originates from the fact that quantitatively accurate and reliable (within 
the accuracy of the approximation for the exchange-correlation part of the total energy of 
the electronic subsystem in density functional theory (DFT) effective cluster interac- 
tions can be obtained only by the Connolly- Williams (CW) or structure inversion method 
on the basis of the total energies of a set of specifically chosen ordered structures 
calculated by the so-called full potential (FP) methods, which have no restrictions on the 
form of the one-electron potential and density. However, the structure inversion methods 
become practically unusable in the case of an inhomogeneous system, not only because 
of the large number of basic structures which must be calculated to extract position- 
(layer-, for instance) dependent interactions, but mainly because of the large size of those 
basic structures (supercells) which are needed to factorize a specific interaction inside 
homogeneous parts of the systems, e.g., inside the layers parallel to a surface. 

In this situation there appears to be only one alternative to the structure inversion 
methods: The so-called generalized perturbation method (GPM) proposed by Ducastelle 
et al. J7||| on the basis of the coherent potential approximation (CPA) |8|-|T0|| and for- 
mulated within tight-binding (TB) theory. Later the GPM was generalized in a straight- 
forward manner [O-n for use in ab initio calculations based either on the Korringa- 
Kohn-Rostoker (KKR) method or the linearized muffin-tin orbitals (LMTO) method in 
the atomic sphere approximation (ASA). The main idea behind the GPM is to calculate 
perturbatively the total energy difference between the alloy in the initial state, which 
is completely random, and in a final state in which only one specifically chosen atomic 
distribution correlation function or short-range order parameter is different from that in 
the random state. This makes the GPM very efficient and convenient to use as it directly 
yields the needed effective cluster interactions. 

However, it is known, although rarely mentioned in the literature, that the interac- 
tions obtained by the GPM yield a quantitatively poor description of the ordering in 
real alloys (see, for instance, Ref. [pj|]) in those cases where there is a substantial size 



mismatch between the alloy components. This failure may only partly be attributed 
to lattice relaxation effects. Rather, it originates not from the GPM method itself but 
is a consequence of inappropriate approximations in the basic methods underlying the 
GPM calculations. This is so, because, as has been demonstrated by Bieber et al. 
in parametrized tight-binding calculations and by Singh et al. JL7| in ab initio KKR-CPA 
calculations, the GPM interactions may provide (under certain conditions) a consistent 
description of the ordering or configurational energy. That is, the ordering energy ob- 
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tained from the GPM interactions, calculated in the framework of a particular technique, 
agrees reasonably well with the ordering energy obtained directly from the total energy 
calculations by the same technique. 

The ab initio techniques underlying GPM calculations are usually the KKR-ASA, 
KKR- ASA-CPA, and LMTO-CPA methods [|18j which are based on a number of approxi- 
mations, such as the CPA, the single-site (SS) approximation for the electrostatic part of 
the DFT problem, and the spherical approximation for the form of the potential which, 
depending on the geometry, is called either the muffin-tin (MT) or atomic sphere ap- 
proximation (ASA). The question is which of these approximations is the most severe in 
the cases where the alloy components have a substantial size difference. To answer the 
question, we note that a size difference leads to so-called "charge transfer effects" or, to 
be more precise, to a non-zero net charge for each alloy component inside their atomic 
spheres chosen to be of equal size. 

Although there are systems, where the CPA may lead to substantial errors, it is clear 
from a general point of view that the CPA cannot be responsible for the errors in the 
case of pronounced charge transfer effects, because the error of the CPA is mainly related 
to specific features of the electronic structure of the individual alloy components such as 
the difference in the position and overlap of the energy bands Moreover, there is a 

number of different calculations which show that in such systems as, for instance, CuPd 
and CuAu, where the alloy components have similar electronic structures but different 



atomic size, the CPA works fairly well p0|-|22|| as a method for obtaining the average 
electronic structure of random alloys. 

As far as the underlying KKR-CPA or LMTO-CPA methods are concerned, much 
larger errors may in fact come from the use of the single-site approximation in the self- 
consistent DFT part of the calculations as this yields no information about the distribution 
of the charge outside the individual atomic spheres of the alloy components. In fact, 
the effective medium outside the individual atomic spheres of the alloy components is 
electroneutral, and therefore, if the net charge of an atomic sphere is non-zero, Poisson's 
equation cannot be solved properly. A number of different models have been proposed 
to include the missing screening charge in the solution to Poisson's equation p3|-p6| and, 
most recently, a general formalism of screened Coulomb interactions (SCI) based on the 
knowledge of the spatial distribution of the screening charge around an impurity has been 
developed in Ref. |27j together with a formalism for the SCI contribution to the GPM 
interactions [0]. Although the SCI may now be included in SS-DFT-CPA calculations, 
this does not solve all problems connected with the description of the energetics of alloys. 

It is not surprising that the main source of inaccuracy in the KKR-CPA and LMTO- 
CPA methods is the spherical approximation, MT or ASA, for the form of the electron 
density and potential (in the following we will consider only the ASA, since the difference 
between the ASA and MT is unimportant for the later discussion and results). For 
instance, in the extreme case, where one of the alloy components is a vacancy, the error 
due to the use of the ASA is about 100% (or several eV in absolute values) for the 
vacancy formation energy [p8|| . As has been shown by Korzhavyi et al. ||2Q] , this kind 
of error originates from the oversimplified description of the non-spherical electrostatic 
contribution to the one-electron potential and energy mainly from the charge density on 
the atoms next to the vacancy. 
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This is similar to the case of surfaces where the quite large ASA error may be substan- 
tially reduced by the inclusion of the multipole moments of the electron charges inside the 
atomic spheres [[30]. The so-called ASA+M approach significantly improves vacancy and 
defect formation energies [^,0, surface energies |32| and alloy energetics [|33|| . Recently, 
Finnis et al. have included the multipole moments in their self-consistent tight-binding 
model which allowed them to obtain quite accurate description of the energetics of zirco- 
nia. 

In this paper we show that the use of the ASA+M approach leads to a representation 
of the configurational part of the total energy, which is very close to the full-potential 
results. Since the polarization of the atoms in an alloy is almost entirely determined 
by their closest local environment, it is obvious that the effect of polarization cannot 
be described properly in the single-site approximation. Nevertheless, the SS-DFT-CPA 
methods may be still used for the electronic structure and total energy calculations of 
random alloys if the definition of the SCI is modified. It is the main purpose of the present 
paper to demonstrate how this may be done in the cases of ordinary bulk homogeneous 
random alloys and inhomogeneous systems such as surfaces. 

The paper is organized as follows. In section II we introduce the ASA+M approxima- 
tion and outline some details of our calculations. In section III we compare the ordering 
energies of NiPt alloys, calculated by the KKR method in different approximations and 



by the Vienna ah initio simulation package (VASP) [p5| , p6| . In section IV we define the 



on-site screening Madelung potential, which should be added to the one-electron poten- 
tial in the SS-DFT calculations. The SCI and the Madelung energy of a random alloy 
are defined in section V. In section VI we calculate the intersite SCI in NiPt in different 
approximations and by different methods. In section VII the screened generalized per- 
turbation method interactions are calculated and compared with the Connolly- Williams 
interactions. In section VII we show how the formalism for the SCI should be modified 
in the case of inhomogeneous systems where there are several non-equivalent sublattices, 
like partially ordered alloys or surfaces. 



II. BEYOND THE ASA 

In a previous publication []27[], in the following referred to as I, we presented a consistent 
and variational, within DFT, approach to the electrostatic screening effects in random 
alloys, and within the ASA we found that these screening effects were almost independent 
of alloy composition, lattice spacing, and crystal structure. However, in those cases where 
the alloy components have a substantial size difference one cannot obtain a quantitatively 
correct description of the configurational part of the total energy of metallic alloys within 
the ASA. One must therefore go beyond the spherical approximation at least for the 
density. In the following we will show how this may be done. All other detail of our 
approach may be found in I. 
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A. Multipole correction to the atomic sphere approximation: ASA+M 



The idea behind the multipole correction is simply to include those contributions to 
the electrostatic multipole-moment expansion of the intercell or Madelung part of the one- 
electron potential and energy, which are neglected in the ASA. If the multipole moments 
of the electron charge, g R , inside an atomic sphere centered at R are defined as 

9r = J s (^~) n R( r n)*L(r R )dr R - Z ^oi, (1) 

where L is short-hand for the (I, m) quantum numbers, Sr the radius of the atomic sphere, 
tir the non-spherical charge density, and Yl a real harmonic, the Madelung contribution 
to the one electron potential is given by 

4^ E M R,R'9r', (2) 

D R',L' 

while the Madelung energy which now includes the multipole-multipole electrostatic in- 
teractions between different lattice sites may be written 

e m = ^e ik E <:r4'- (3) 

ZD R,L R',L' 

L L' 

In these expressions, M R ' R , is the multipole Madelung matrix which is equivalent to the 
conventional LMTO structure constants and the number of multipoles included in the 
L, V summations is determined by the angular momentum cutoff l max in the basis set used 
in the Green's functions calculations. Owing to the properties of the Gaunt coefficients 
non-zero multipole moments of the charge density may be generated for /-values up to 

We note that, since in the ASA+M the one-electron potential is still kept spherically 
symmetric inside each atomic sphere, the only term which contributes to the one-electron 
potential is the L=(0,0) or monopole term. This simple restriction on the form of the-one 
electron potential violates the variational connection between the Madelung potential and 
energy and, in turn, between the one-electron potential and the total energy, i.e., 

v M r = % R ^ -? — ■ 4 

However, since this is just a consequence of the model, but not of theory in general, it 
does not create any problems. On the other hand, the reinstatement of the variational 
connection between the one-electron potential and the total energy by keeping only the 
monopole-multipole term in (||) may lead to a substantial error in the total energy, as the 
multipole-multipole interactions will not be accounted for. 



B. Details of calculations 

The Green's function technique has been used in both the KKR-ASA and the locally 
self-consistent Green's function (LSGF) calculations in the scalar relativistic and atomic 
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sphere approximations. This part of the techniques is described in Refs. p7| , |32|j37| . The 
basis functions have been expanded up to l max = 3 (spdf-basis) inside the atomic spheres, 
while the multipole moments have been calculated up to l^ ax = 6. We have also performed 
a number of calculations in the spd-b&sis, i.e., l max = 2 and l^ ax = 4. The integration 
of the Green's function over energy was performed in the complex plane over 16 energy 
points on a semicircular contour using a Gaussian technique. We have used the generalized 
gradient approximation (GGA) of Perdew and co-workers [B8J. For each structure the 
integration over the Brillouin zone (BZ) has been done by using equally spaced k-points 
in the irreducible part of the appropriate BZ and the number of k-point has been chosen 
to be equivalent to 500 - 1000 uniformly distributed k-points in the irreducible part of the 
BZ of the fee structure. Core states were recalculated at each DFT iteration. 

For benchmark calculations we applied the Vienna ab initio simulation package 
(VASP) described in detail in Refs. p5|p6| . These calculations were performed in a plane- 
wave basis, utilizing fully non-local Vanderbilt-type ultrasoft pseudopotentials (US-PP) 
| 3~9fl which allow the use of a moderate cutoff in the construction of the plane- wave basis 
for the transition metals. In the actual calculations the energy cutoff was set to 302 eV, 
exchange and correlation were treated in the framework of the GGA and the inte- 
gration over the Brillouin zone was performed on a Monkhorst-Pack k-mesh [ |4U[| . Test 
calculations showed that, depending on structure, the required convergence was reached 
for 35 to 275 k-points in the irreducible wedge of the BZ. 



III. ORDERING ENERGIES IN NI-PT 

We start by demonstrating the accuracy of the various approximations which are 
usually used in KKR(LMTO)-ASA-like calculations. For this purpose we have calculated 
a set of ordered fee NiPt alloys by the KKR-ASA method and by the US-PP for a fixed 
lattice constant without any local or anisotropic relaxations. The lattice constant has 
been chosen to be a ~ 3.791 A, which corresponds to an atomic Wigner-Seitz radius of 
2.8 a.u.. The ordered structures include: Ll 2 (Cu 3 Au-type), D0 2 2 (TiAl 3 -type), Z3, (3 
4TJ, 7 (Pt 2 Mo-type), Ll (CuAu-type), CH or "40" (NbP-type), Z2 |4j, Ll x (CuPt-type), 



and the so-called SQS-16 

To simplify the comparison we present in Table [I] the values of the calculated " mixing" 
energies of the above mentioned ordered structures 

^ ZpNi _i_ „ TpPt 

j?Ni m Pt n _ j?Ni m Pt n _ ,nCj tot ' nCj tot (c\ 

m + n 

where E^l and Ef* t are the total energies of the pure components calculated at the same 
lattice constant. All total energies are per atom. The mixing energies of the random 
alloys have been obtained on the basis of all the energies included in Table |, except 
SQS-16, plus the energies of Ni 7 Pt and PtyNi (CuPty-type) ordered alloys (not presented 
in the Table) by the Connolly- Williams method in which the total energy expansion 
included pair interactions at the first four, 7th and 10th coordination shells (these are the 
largest pair interactions in this as well as in many other fee transition metal alloys), the 
first four triangle interactions and the two tetrahedron interactions corresponding to the 
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tetrahedra of nearest neighbors, and the straight line along [111] direction (the last being 
quite substantial in many systems). 



The SQS-16 is a so-called special quasirandom structure [M\ which consists of 8 atoms 
of one type and 8 atoms of another type distributed in the unit cell in such a way that 
the first 7 pair, the nearest neighbor triangle, and the tetrahedron atomic distribution 
correlation functions are the same as in the random alloy. Hence, the fact that the values 
of E^x~ m an d EmS? are nearly eaqual indicates that: 1) the SQS-16 provides a good 
model for the random NiPt alloy, and 2) the convergence of the CW method is reasonably 
good. Part of the convergence of the CW method is provided by the use of total energies 
on a fixed lattice, whereby the volume dependent contribution to the total energy is not 
expanded in terms of cluster interactions, which is an ill-defined procedure in metallic 
systems and usually leads to very bad convergence of the CW method (see, for instance, 



In Table [I] we present KKR-ASA results in the spdf- as well as the spd-hasis and in 
both cases we show results in the pure ASA, i.e., without multipole moment contributions, 
and in the ASA+M. We also include the results of neutral sphere calculations, in which 
the atomic sphere of Pt is chosen to be larger than that of Ni in order to provide zero net 
charges of the atomic spheres. Although the comparison of the mixing energies should be 
done with some caution, because the ground state properties of the alloys are different 
in different approximations and because all the calculations have been performed at the 
same fixed lattice constant, it is clear that the ASA+M approach in the spdf -basis leads 
to values of E m i X which are in considerably better agreement with the US-PP results than 
any of the other approaches. 

Using the results of the Connolly- Williams method for the total energy of random 
alloys one can calculate the ordering energies, defined as the difference between the total 
energies of an ordered and a random alloy at the same composition. In Table |I| we 
compare the ordering energies of different structures calculated in different approximations 
for NiaPt, NiPt, and PtaNi. Such a comparison makes sense since the ordering energies 
are much less volume-dependent than the total energies themselves. Again, it is seen that, 
relative to the US-PP results, the ASA+M approach in the spdj '-basis gives not only the 
best, but also a quite accurate description of the ordering effects in NiPt. It is also seen 
that the ordering energies in the ASA+M but without multipole-multipole contribution 
to the Madelung energy appear to be half way beteeen ASA and ASA+M results. As 
we will see below, this is in fact the limit of accuracy which can be reached in consistent 
single-site mean-field calculations. 

It is also obvious from the table that the KKR-ASA does not yield reasonable values 
for the ordering energies in the case of transition metal alloys unless /-states are included 
in the basis. These /-states are needed to supply a better augmentation of the basis 
functions at the atomic sphere and a better interstitial charge density P5| . The neglect 



of /-states can only be partly compensated by the use of the so-called combine-correction 
term in the LMTO method . 



Another important conclusion, which can be drawn from the results in Table [II], is 
the fact that the use of neutral spheres leads to a substantial underestimate of the order- 
ing effects in KKR-ASA(+M) calculations. In other words, although the neutral-sphere 
approach formally solves the problem of the electrostatic interaction in an alloy, the elec- 



7 



trostatic contribution to the one-electron potential and energy being zero by definition, 
it introduces errors which are unacceptable in a quantitative description of the configura- 
tional energetics. 

The reason for this failure is the following. If we compare the values of E mix from 
Table | obtained with neutral spheres with those obtained with equal spheres, we find 
that the neutral sphere approach leads to substantial lowering of the total energy of the 
ordered alloys. However, the amount of lowering is structure dependent: In the case 
of "phase separated" -like structures, where there is a certain clustering of atoms of one 
type, like Z3, Z2, the lowering is much greater, than it is in some of the more "normal" 
structures. Through the CW procedure the exaggerated lowering of the energies of the Z3 
and Z2 structures leads to an exaggerated lowering of the total energy of the random alloy, 
making, in the end, all the ordering energies much smaller in absolute value, than they 
should be. In the case of the "phase-separated" -like structures this kind of error originates 
in the poor filling of the crystal space provided by the non-equal atomic spheres. The 
division of a crystal into atomic spheres is in itself quite a crude approximation, the error 
of which is not generally known. Now, if such division is done differently for different 
ordered structures by means of non-equal atomic spheres, it is obvious that the resulting 
errors will be structure dependent and render a comparison of total energies meaningless. 

Based on the above considerations it appears that the use of equal atomic spheres for 
the alloy components provides the only consistent and correct way of dividing space for 
configurational or ordering energy calculations, if the underlying lattice is a simple Bravais 
lattice, such as fee and bcc. In that case the packing of space by the atomic spheres of the 
alloy components is homogeneous and independent of the alloy configuration which makes 
at least part of the ASA error systematically compensated when the energy difference of 
different structures is calculated. Of course, the situation is different in alloys where 
the underlying lattice has a more complicated structure. However, in the configurational 
energy calculations one should follow the recipe of choosing equal sphere radii of the 
alloy components on the sublattice, where the alloying (ordering) is taking place (see, for 
instance, p7| ). 

IV. THE ON-SITE SCREENING MADELUNG POTENTIAL IN THE ASA+M 

Although the ASA+M approach seems a natural generalization of the ASA, that 
changes little in the formalism, it has a large effect on the way the SCI must be treated 
in random alloys. The most dramatic consequence of the ASA+M is the fact that the 
correct SCI can no more be obtained in single-site electronic structure or Green's function 
calculations: The largest polarization effects, which give multipole-moment contributions 
to the SCI, actually comes from the nearest neighbors of the atom the screening of which 
is considered. 

This is easy to see in the case of a single impurity in an otherwise perfect crystal which 
is the dilute limit of a random alloy. The point-group symmetry of the impurity site is 
exactly the same as that of the underlying lattice, while none of its neighbors has even 
inversion symmetry and, therefore, every atomic sphere around the impurity has a non- 
zero dipole moment. Such a dipole moment plus the higher multipole moments induced 
by the presence of the impurity can only be found in Green's function calculations which 
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include these neighboring sites in the perturbation part of the Dyson equation. Thus, the 
multipole-moment contribution to the SCI can be obtained only beyond the single-site 
approximation in the Green's function calculations. 

Nonetheless, we will determine a simple parameterization of the SCI for the SS-DFT- 
CPA method, which will allow us to obtain accurate results in SS-KKR(LMTO)-CPA 
calculations. Such a parameterization of the SCI in the ASA+M is purely a fitting pro- 
cedure in contrast to the SCI determined by SS-ASA-Green's function calculations in I 
basically for the purpose of showing the existence of a consistent theory of the SCI in 
random alloys in the framework of the single-site- CPA theory. 

Let us first consider the on-site Madelung potential in a random alloy due to the 
screened Coulomb interactions which should be used in SS-DFT-CPA calculations for the 
random alloy. In the case of a binary random alloy, Ai_ c B c , this potential can be defined 
as (see I) 

Vi = -e 2 a scr (0)^, (6) 

where qi is the net charge of the zth alloy component, S the Wigner-Seitz radius, and 
a l scr (0) the on-site screening constant which can be obtained from the screening charge by 
performing "impurity" calculations (see I), for instance, by exchanging the type of atom 
at a particular site of the supercell modeling the random alloy. Then a scr (0) is 

2S 

cw(0) = — E S ,l(R)Ql(R), (7) 

e L,R 

where Si,jy(R/) are the canonical structure constants, the first few terms of which are 



explicitly defined, for instance, in Ref. | 48fl . Further, Ql(R) are the normalized multipole 

ing che 

Ql(R) 



moments of the screening charge in the atomic sphere centered at R 

Aq L (R) _ Aq L {K) 



Ag(R=0) Ag L=00 (R = 0)' 

where Ag^(R) is the difference between the L- moments of the charge in the atomic sphere 
at R) after and before the impurity has been introduced at the site R =0. 

There is one important point. In the pure ASA a scr (0) can be determined in an 
alternative although formally equivalent manner by 

S ^ D' ^ 

aw(0) = (a rand ) = — — , (9) 

e 2 < q { > 

where < Vi > and < > are the average values of the Madelung potential and net charges 
of the zth alloy component in the self-consistent supercell calculations of the random alloy. 
Being practically exact in the ASA this scheme is only approximately valid in the ASA+M, 
where the average Madelung potential of the supercell, which is equal to the Madelung 
potential of the underlying lattice, is not equal to zero in general due to the presence of 
the non-zero multipole moments at least for I =4. In the case of inhomogeneous systems 
like surfaces, where the multipole contribution to the average Madelung potential is quite 
large, Eq. (||) cannot be used at all. Therefore, @ should be modified by substracting 
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the corresponding average values of the net charges and Madelung potentials of the alloy 
components on the (sub)lattice, q = (1 — c) < qA > +c < qs > and v = (1 — c) < va > 

+c < v B > 

, N S < v i > —v . . 

a SC r(0) = (ot rand ) = — -. (10) 

e 2 < qi > -q 

In Fig. [I] we compare the monopole moment of the normalized screening charge, 
Ql=oo(R), in a Ni 50 Pt5o random alloy obtained by changing the type of the atom from Pt 
to Ni at a some particular site in a 384-atom NisoPtso supercell, the first seven short-range 
order (SRO) parameters of which are equal to zero. The calculations have been carried out 
in the single-site (SS)-LSGF-ASA as well as the embedded-cluster (EC)-LSGF-ASA+M 
methods (see I). The local interaction zone (LIZ) in the EC-LSGF-ASA+M calculations 
included 43 atoms, i.e., the central atom and its first three coordination shells (LIZ=4), 
while in the SS-LSGF calculations the LIZ included only one atom (LIZ=1). 

It may be seen from Fig. [I] that the screening is more efficient in the ASA+M than 
in the ASA, although the difference between the two cases is very small. Nevertheless, 
the effect of the multipole moments on the on-site screening constant a scr (0) is quite pro- 
nounced, it increases to about 0.74 from the 0.61 in the ASA calculations. The largest 
multipole contribution comes from the dipoles in the first coordination shell, which con- 
tribute almost 0.1 to the a scr (0), while the quadrupole and octuple moments contribute 
0.036 and 0.016, respectively. 

The effective charge transfer, Aq =< q^i > — < qp t >, increases from 0.505 in the 
ASA to 0.583 in the ASA+M. However, it is still reproduced correctly in SS-DFT-CPA 
calculations, provided the correct, i.e., corresponding to the ASA+M, value of the on- 
site screening constant, a scr (0), is used in (^). This is so, because the on-site screening 
constant is a parameter which determines a constant shift of the one-electron potential 
both in the ASA and ASA+M. This constant may be chosen to contain the monopole- 
multipole interactions of the charge at a given site with its screening cloud, thereby 
yielding the correct effective charge transfer. Unfortunately, the same is not the case for 
the Madelung energy and the SCI in general. 



V. THE SCI AND MADELUNG ENERGY OF A RANDOM ALLOY 

The screened Coulomb interactions, V scr (R), are the energies of the electrostatic in- 
teraction between the electron charge density inside an atomic sphere centered at some 
site (it is convenient to choose this site as the origin) and the perturbed electron density 
and its screening charge at another, in general, different site. In the ASA+M, the SCI 
can be represented in the form of the multipole expansion, 

v scr (R) = Y,Vs L cr(R) , (ii) 

L 

where U s ^ r (R) is the L-component contribution to the SCI at the distance R, which for a 
binary A!_ C B C alloy on a Bravais lattice may be expressed as in I 
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''-Aq 2QLiK 



o)«L(R) 



e 

45 



Ag 2 Q L (R 



5* 

= 0) X! ^l=oo,l'(R')Ql'(R'-R) 



;i2) 



Here Ag is the difference of the net charges in the atomic sphere after and before the 
perturbation, i.e., the exchange of the type of atom at site R, af cr (R) is a generalized 
screening constant, Sl,l'(R') are the canonical structure constants |4"8|| , Qx(R'— R) are the 



normalized multipole moments of the screening charge, and AqL =00 =< qA > — < qB > 
is the effective charge transfer in the alloy. 

Within the single-site mean-field considerations presented in I, all the multipole mo- 
ments on the alloy sites are uncorrelated, the average value of Ql(R = 0) being either 
very small or equal to zero, unless L = 00, and thus the only non-zero SCI is V s ^r 00 , 
which, for instance, in the case of the fee underlying lattices can be written in the form 
(see I) 

V scr (R) = A ? 2 ^^P, (13) 



2 S 



where 



a SC r{R) = \ £ S l= oo,l>(R')Qv(R , -'&). (14) 

The on-site term V scr (R = 0) is the energy of the electrostatic interaction between the 
net charge of an alloy component and its screening density or, as has been shown in I, it 
is the screening Madelung energy of the random alloy. It is easy to see that this energy 
is DFT-consistent with the corresponding screening Madelung shift of the one-electron 
potential @. However, in contrast to the screening Madelung potential which correctly 
reproduces the effective charge transfer in random alloys, V scr (R = 0) underestimates the 
corresponding Madelung energy in the ASA+M because of the missing contribution from 
the multipole-multipole interactions in the single-site mean-field approximation. In the 
case of, for instance, a Ni 50 Pt5o random alloy this contribution is about —3 mRy/atom 
for the Wigner Seitz radius S =2.8 a.u.. 

This means that if one wants a quantitatively accurate value of the total energy of 
a random alloy in the SS-DFT-CPA calculations consistent with the supercell ASA+M 
calculations one needs to modify the definition of the Madelung energy of the random 
alloy by introducing a fitting parameter. The simplest way to do so is to define the 
Madelung energy of the random alloy as 

^ = E c ^r, (is) 

i 

where c, is the concentration of the zth alloy component, and 

Er = jp q f scr{R s =0 \ (i6) 
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which means that in the case of a binary random alloy 

E r ^ n a d d =c(l-c)(3V scr (R = 0). (17) 

Here, f3 is the renormalization coefficient which is approximately equal to 1.16 for most fee 
and hep transition metal random alloys. Thus, E? cr and V % scr are no more DFT-consistent 
(see®). 

As discussed above, this violation of general theory is a consequence of the ASA+M, 
which on the other hand brings the ordering energies of the much more efficient SS-DFT- 
CPA approach into good quantitative agreement with the corresponding full-potential 
results. Although the difference between the multipole-multipole and monopole-multipole 
results in Table [TI] might not look so dramatic, the omission of the multipole-multipole 
interactions in the ASA+M Madelung energy has much more serious consequences in the 
case of, for instance, surface energy anisotropy calculations, which cannot be reproduced 
even qualitatively without this term. 

Finally, we show that the CPA itself introduces relatively small errors in the elec- 
tronic structure calculations of a random alloy. In Fig. we compare the total density of 
states (DOS) of a random Ni 50 Pt5o alloy and the local, Pt and Ni, contributions obtained 
by two different methods, SS-KKR-CPA and EC-LSGF. The SS-KKR-CPA calculations 
have been performed with a scr = 0.74 for the screening Madelung potentials of the alloy 
components @, and in the EC-LSGF calculations LIZ=4 has been used to calculated 
the SQS-16 supercell of Ni 50 Pt 5 o. Since, all the SRO parameter of the SQS-16 are equal 
to zero up to the 7th coordination shell, beyond which the net charge at each site is 
practically completely screened, the EC-LSGF results may be considered as a benchmark 
for the SS-KKR-CPA method. The agreement between the EC-LSGF and SS-KKR-CPA 
results is seen to be good clearly indicating that the CPA works fairly well for this system. 



VI. INTERSITE SCREENING CONSTANTS IN NIPT 

There are several ways of obtaining the intersite screening constants, a scr (R), which 
determine the corresponding intersite SCI. First, they may be obtained directly from the 
normalized moments of the screening charge, Ql(R), by menas of (|14"1). This requires two 
self-consistent supercell calculations: One performed for some initial atomic configuration 
in the supercell and the other for the same supercell with an "impurity" at the site where 
the type of atom is changed (see paper I). In this manner, however, only the monopole- 
multipole part of the SCI can be found. 

Secondly, one may take advantage of the special properties of the LSGF method, 
namely of the fact that the electronic structure in the LSGF method is obtained in the 



so-called combined cluster-effective medium approach |27],|37]]. That is, the local environ- 
ment effects are taken into account only inside the LIZ during the electronic structure 
calculations, while the rest of the crystal is seen by each atom as a random alloy, described 
by the CPA effective medium. This means that all the correlated atomic configurations 
attributed to the non-zero SRO parameters beyond the LIZ do not contribute to the elec- 
tronic structure of the supercell. At the same time, since the Madelung problem is solved 
exactly, they are accounted for in the electrostatic part of the total energy. 
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Considering the change of the atomic configuration on the lattice as a small pertur- 
bation, one may argue on the basis of Andersen's force theorem [4j| that the difference 



in the total energies of the alloy with some non-zero SRO parameters and a completely 
random alloy in the LSGF calculations may be given by 



rpLIZ=N — jpSRO _ jprand 
^ord — ^tot 111 tot 

1 



2V-1 N scr 

aiZiV one - e i(Ri) + Y aiZiV scr (Ri 

i=l i=l 



where ctj are the Warren-Cawley SRO parameters (see, for instance, Ref. 0) at the zth 
coordination shell, N the size of the local interaction zone, which is determined as the 
number of coordination shells around the impurity site plus 1, and N scr the coordination 
shell beyond which the SCI vanish. In the case of an fee alloy it may safely be assumed 
that, N scr = 7. It is obvious that the first term in flT8| ) may be identified with the usual 
GPM interactions. 

It follows from (|T%|), that if iV — 1 < N scr and «j=0 for all i < N scr except for one 
coordination shell, j, which is beyond the LIZ, j > N — 1, then the intersite SCI for this 
specific coordination shell can be determined in two LSGF calculations: One performed 
for a supercell which corresponds to the random alloy (ctj = for all i < N scr ) and the 
other for a supercell which satisfies the above described conditions with ctj ^ 0. In this 
case 

= e*c(l-c)°A d q ^*i (19) 
In the ASA all the SCI can be determined by using (|T§[) in the corresponding SS- 



LSGF-CPA calculations, while a scr (R2) is the first SCI which can be determined in the 
ASA+M in this way, since there is no multipole contribution to the SCI in the single- 
site approximation (N(LIZ) = 1). It is also clear that if the multipole-multipole Madelung 
energy is included in the corresponding LSGF calculations, then, in principle, a scr (R) 
determined from ( |19| ) will also contain the multipole-multipole contribution. 

However, in this approach the higher order atomic distribution correlation functions 
(multisite SRO parameters) of the supercell should be also optimized. This is so, since all 
the multisite interactions for the figures inscribed in the LIZ with an vertex located at the 
central atom of the LIZ, also contribute to the ordering energy. Moreover, in the case of 
the ASA+M, in principle, there is a non-zero contribution from the SCI to the multisite 
interactions themselves. In our calculations, we have not optimized the multisite SRO 
parameters of the supercells and, therefore, considering the quite small values of the SCI 
beyond the first coordination shell we have not used (|19"D in the ASA+M calculations. 

In Table [TTI| we compare a scr (R) for the first four coordination shells calculated either 
by direct summation of the normalized multipole moments fllil) or by means of (|19D in 
the SS-LSGF total energies. The LSGF calculations have been performed for 384-atom 
supercells of an equiatomic Ni 50 Pt 5 o alloy with the corresponding sets of SRO parameters. 
Although we have not optimized the higher order SRO parameters, they turned out to be 
small: ~ 0.02 or less, at least for the triangle and tetrahedra of the nearest neighbors on 
the fee lattice. Thus they should not affect the results, at least in the ASA. As seen from 
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Table [TTJthe agreement between the intersite screening constants, a scr (R), determined in 
two different calculations in the ASA, is quite reasonable. 

As shown in I, a scr (R) is practically a universal function for metallic alloys on simple 
Bravais lattices in the ASA and the single-site approximation. In the ASA+M approach, 
this is not the case any more. However, our calculations for the dilute limit of NiPt alloys 
and for different lattice constants show that a scr (R) changes very little in these cases and, 
in fact, less than the difference between the values of a scr (R) calculated in two different 
ways. For instance, in Table |TTT| we show the values of a scr (R) obtained for a Ni impurity 
in pure Pt at S — 3 a.u., which are very close to those of the random equiatomic alloy at 
S = 2.8 a.u. 



VII. SCREENED GPM INTERACTIONS IN NIPT 

As demonstrated in I the intersite screened Coulomb interactions, V scr (R), given by 
(IHf) must be added to the corresponding one-electron term given by the GPM in order 
to satisfy the force theorem In the case of a binary random alloy A c Bi_ c the 



screened generalized perturbation method (SGPM) interactions are defined as in Ref. 

V SGPM (R) = V GPM (R) + V scr (R) (20) 
= V GPM (R) + ^ Aq 2^S. : 

where V GPM (R) is the usual GPM interactions obtained as the change of the one-electron 
energies due to specifically induced alloy configurations on the alloy's underlying lattice. 
In fact, V GPM (R) should be renormalized due to the intersite SCI. However, this problem, 



as well as the complete SGPM formalism will be considered elsewhere |5T[ . 

In Table [IV] we compare the first four most important effective pair interactions (the 
rest of the pair interactions and the multisite interactions are less than 0.1 mRy) in the 
NiPt fee alloys at S =2.8 a.u. obtained by three different techniques, i) the SGPM for 
equiatomic alloy composition, ii) the Connolly- Williams method on the basis of the total 
energies of the ordered alloys described above, and iii) direct calculation from the EC- 
LSGF total energies of Ni 50 Pt5o alloys similar to the case of the intersite screening constant 
calculation described in the previous section. That is, in iii) the effective interactions 
have been obtained from the EC-LSGF total energy calculations for a completely random 
alloy, cti =0 for all i < N scr , and for an alloy with one non-zero aj. However, now 
j < N — 1, where N is the size of the LIZ in the LSGF calculations, and therefore 
the local environment effects attributed to the non-zero SRO parameter, are included in 
the electronic structure calculations. In this case the effective pair interactions can be 
determined as 

o rpLIZ=N 

V 3 = — . (21) 

C(l — C)ZjO!j 

It is important, that (|2l| ) is not based on any additional approximations and therefore 
constitutes a direct way of determining the effective interactions with an accuracy which, 
in principle, is restricted only by the approximations used in the LSGF calculations, that 
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is, mainly by the ASA+M, since the CPA yields very small relative errors in the EC-LSGF 
calculations with LIZ> 2. 

Although the SGPM interactions as well as the interactions determined from ( pT| ) are 
concentration dependent, while the Connolly- Williams interactions are concentration in- 
dependent, they can be compared since the concentration dependent interactions obtained 
for an equiatomic alloy composition are equal to those of the concentration independent in- 
teractions, at the same fixed volume p3| . Therefore, if the basis in the Connolly- Williams 
method includes all the important interactions for a given system, the Connolly- Williams 
interactions obtained from the KKR- ASA+M calculations of the ordered alloys at a fixed 
lattice constant should be equal to those obtained from ( pip for the equiatomic alloy 
composition at the same lattice constant. 

It is clear from Table IV], that in general the agreement between the Connolly- Williams 
and the direct calculations is quite good. Further, the Connolly- Williams interactions 
obtained by the KKR-ASA+M and by the US-PP (CW-KKR-ASA+M and CW-US- 
PP) agree well with each other, except for V3 which is a little larger in the CW-US-PP 
calculations than in the CW-KKR calculations. This, then, confirms our point that the 
ASA+M approach allows us to obtain quantitatively accurate configurational energies of 
metallic alloys. 

It is interesting to note that although the neutral-sphere approach yields quite large 
errors for the ordering energies, it seems to work remarkably well for the GPM interac- 
tions, except for V±, which is more than twice as small as it should be. It is probably a 
coincidence that they come out very close to the CW-KKR-ASA+M interactions, since 
the SCI are quite small beyond the first coordination shell. Unfortunately, the CW in- 
teractions obtained in the neutral sphere calculations do not seem to be convergent, the 
three- and four-site interactions being of the same order of magnitude as the interaction 
at the first coordination shell and therefore, they are not given in the table. 

The SGPM interactions, obtained in the ASA+M and in the ASA (in these two cases 
we have used the intersite SCI calculated by fllTf)) are not very different, except for the 
interactions at the first coordination shell. This is most probably due to the missing 
multipole-multipole contribution to the intersite SCI in the ASA+M. One may also see 
from Table Table |Vj, that the SGPM-KKR-ASA+M interactions are in fact quite close to 
those of the CW-KKR-ASA+M obtained without multipole-multipole electrostatic inter- 
actions (0-L). These interactions have recently been used in Monte Carlo simulations of 
the ordering in NiPt and reproduced quite well the order-disorder transition temperature 
for an equiatomic alloy composition and the values of the SRO parameters in a random 
alloy at T=1200 K fSJ. 

Finally, in Table [V] we show the ordering energies of the four equiatomic ordered 
alloys, obtained from the first 20 SGPM interactions, although the contribution from the 
interactions beyond the fourth coordination shell is only a few percent of the total ordering 
energy. Comparing these energies with those from the direct total energy calculations, 
presented in Table [TI], we find reasonable agreement for the ASA+M and the ASA results 
and very good agreement in neutral-sphere approach. 
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VIII. THE SCI AT ALLOY SURFACES 



The generalization of the SCI formalism to inhomogeneous systems, such as partially 
ordered alloys or surfaces, is straightforward. In the latter case the SCI also become 
inhomogeneous and, therefore, definition (|T3|) should be rewritten as 



V^(R) = jA qx A qy ^^, (22) 

where Aq x is the effective charge transfer on the A sublattice, and a^(R) the screening 
constant, which is defined as 

a£(R) = l E S^otwCRO^R'-R)- (23) 

Here, the vector R connects the site on sublattice A, where the perturbation of the charge 
density is induced to the site on sublattice A', at which the SCI is determined (see paper I 
for details), and Q^(R) are the multipole moments of the screening charge in the atomic 
sphere centered at R normalized by Aq x (R = 0). 

It is clear that a^(R) depends on the direction of R, i.e., a^(R) ^ a^{R), and 
therefore the SGPM interactions, which in the inhomogeneous systems should be invariant 
under a sublattice index interchange, are defined as 

V x f PM (R) = V^ M {R) + \{V^{R) + Vl\R)) = (24) 

= V^ M {R) + jAq x Aqy a ^ R) + a ^( R \ 

where V x G x f M (R) is the usual GPM interaction. 

As an example of the inhomogeneous system we have chosen three low-index fee sur- 
faces: (111), (100) and (110). The calculations has been performed by the EC-LSGF 
method with an inhomogeneous effective medium, which was fcc(lll), fcc(100), and 
fcc(llO) slabs consisting of 3 vacuum and 6 atomic layers, 4 vacuum and 6 atomic lay- 
ers, and 6 vacuum and 10 atomic layers, respectively, parallel to the surface. The actual 
supercells for the impurity calculations were built on the basis of the effective medium 
supercells by N x x N y translations in the plane parallel to the surface, which were 6x6 
for the (111) and (100) surfaces, and 6x4 for the (110) surface. In the case of the (110) 
surface N X =Q was the period in the closed-packed [110] direction, while N y =4 was the 
period in the [001] direction. Since the screening is insensitive to the alloy composition 
and the lattice parameter, the screening density was obtained for a Ni impurity in pure 
Pt at S =3 a.u.. 

In Table [VT| we present the on-site screening constant, a^ cr (0) = a^ r (R = 0), in the 
first three layers of these surfaces (the impurity is in the A layer). It is clear that the 
surface makes the screening more efficient, since the larger value of the on-site screening 
constant means a closer position of the screening charge to the impurity, and, in fact, 
a scr 1 (0) increases from the most close-packed (111) to the most open (110) surfaces. One 
can also see that the surface influences the screening mainly in the first layer for the 
fcc(lll) and (100) surfaces, and very little in the second layer of the fcc(110) surface. 
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The intersite screening constants , a^(R), are presented in Table |VII[ Again, one 



can see that the surface has a quite substantial effect on the screening constants a\^ r {R), 
which is due to the perturbed electron density in the first layer, A =1. However, the 
values for a^ r (R) are already very close to the corresponding bulk values (see Table fTJ) 
for A =2. 



IX. CONCLUSION 

The polarization of the electron density of the alloy components due to their size 
mismatch makes a substantial contribution to the electrostatic energy of the alloy. This 
contribution, which is missing in the pure ASA, may be accounted for in the ASA+M 
approach through the multipole-moment interactions in the Madelung part of the elec- 
trostatic problem, and, as we have shown, it plays a crucial role in obtaining the correct 
ordering energetics. 

We have also demonstrated that the neutral sphere approach based on the use of 
different atomic spheres for the alloy components on the corresponding underlying lattice 
(or subllatice, in general) leads to unacceptable, quantitative errors, and therefore the 
only consistent way of obtaining correct configurational energetics is to use spheres of 
equal radii for the alloy components on the sublattice where the alloying is taking place. 

Since the multipole moments due to polarization effects originate from the specific 
local atomic configuration around each site, they may in principle be accounted for only 
by methods which go beyond the single-site approximation in the electronic structure 
(Green's function) calculations. However, the use of a simple parametrized form for the 
on-site Madelung potential and energy in the SS-DFT-CPA calculations still allows one to 
obtain a reasonably accurate description of the electronic structure (if the CPA works for 
a given system) and total energy, although, obviously, such a parameterization is possible 
only on the basis of the calculations by more accurate methods. 

The monopole-multipole intersite SCI have been obtained for NiPt fee bulk and sur- 
face alloys. The SGPM interactions, which are the usual GPM interactions plus the 
SCI, reproduce the corresponding monopole-multipole KKR-ASA+M results, which give 
a semi quantitatively correct description of the ordering in NiPt. 
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FIG. 1. The normalized screening charge in fee NiPt, obtained in the ASA and ASA+M 
approximations . 




FIG. 2. The total and site-projected density of states in NisoPtso obtained in the supercell 
EC-LSGF calculations and by the SS-KKR-CPA method. 
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TABLES 



TABLE I. "Mixing" energies, — E m i x (in niRy/atom) of ordered and random NiPt alloys 
obtained by different methods at a fixed Wigner-Seitz radius of S=2.8 a.u. 
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22.5 


P 


16.4 


10.7 


27.7 


23.7 


22.5 


19.4 










Ni 50 Pt 5 o 






Ll 


27.5 


23.8 


33.1 


39.0 


38.4 


28.0 


Z2 


13.4 


5.1 


25.1 


17.0 


13.4 


15.0 


"CH" 


26.3 


22.8 


^1 n 

O-L.U 


n 


37.4 


26.9 


Lli 


22.7 


17.5 


32.8 


30.8 


30.1 


23.0 


SQS-16 


20.7 










21.7 


random 


20.6 


14.9 


29.3 


28.7 


26.8 


21.6 










Ni 66 Pt 3 3 






P 


17.6 


10.0 


23.9 


23.9 


20.6 


19.3 


"Pt 2 Mo" 


21.9 


18.1 


26.2 


31.6 


30.2 


22.9 










Ni 25 Pt 75 






Z3 


14.0 


8.4 


18.4 


18.5 


15.8 


14.6 


D0 22 


19.7 


16.7 


21.8 


28.3 


27.0 


20.2 


Ll 2 


20.0 


16.8 


22.0 


28.2 


27.2 


20.6 


random 


15.4 


10.7 


20.4 


21.1 


19.1 


16.4 
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TABLE II. Ordering energies (in mRy/atom) NiPt alloys at a fixed Wigner-Seitz radius of 
S=2.8 a.u.. Ordering energies obtained without multipole-multipole contribution in the ASA+M 
calculations are given in parentheses. 









KKR-ASA 








Alloy 




spdf 




spd 




US-PP 




ASA+M 


ASA 


neut. 


ASA+M 


ASA 










Pt 3 Ni 








Ll 2 


-6.13 (-7.07) 


-7.84 


-4.69 


-9.17 


-10.58 


-5.98 


D0 22 


-5.06 (-5.91) 


-6.63 


-2.67 


-8.34 


-9.04 


-4.63 


Z3 


2.34 (2.72) 


3.06 


0.31 

NiPt 


-3.85 


4.09 


1.83 


Ll 


-6.91 (-8.05) 


-8.92 


-3.82 


-10.31 


-11.63 


-6.49 


79 


i . lo ( o.OO J 


y.oo 


A 1 A 


1 1 79 


lo.oO 


0.00 


"CH" 


-5.66 (-6.86) 


-7.87 


-1.67 


-10.05 


-10.68 


-5.30 


Lli 


-2.03 (-2.38) 


-2.58 


-3.54 


-2.15 


-3.39 


-1.47 


SQS-16 


-0.07 (-0.12) 




Ni 3 Pt 






-0.18 


Ll 2 


-4.58 (-5.48) 


-6.12 


-1.66 


-7.05 


-8.05 


-4.52 


D0 22 


-4.29 (-5.25) 


-6.01 


-1.37 


-7.12 


-7.09 


-4.20 


Z3 


1.35 (1.92) 


2.31 


1.99 


2.65 


3.33 


1.48 



TABLE III. a scr (Ri) at the first four coordination shells in the ASA and ASA+M obtained 
by different methods. 



Method 


i=l 


2 




3 


4 








ASA 






Eq. (13) 


0.1584 


-0.0017 




-0.0163 


-0.0108 


Eq. (18) 


0.1640 


-0.0026 




-0.0189 


-0.0116 








ASA+M 






Eq. (13) 


0.1279 


-0.0023 




-0.0101 


-0.0050 


Eq. (13) (impurity) 


0.1304 


-0.0035 




-0.0106 


-0.0052 
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TABLE IV. Effective interatomic interactions at the first four coordination shells obtained 
by different methods (in mRy). The SCI contribution to the SPGM interactions is given in 
parentheses. 



Approximation 


Vi 


v 2 


v 3 


v A 






SGPM 










ASA+M (0-L) 


14.05 (15.44) 


0.32 (-0.10) 


-1.09 (-1.22) 


-1.76 


(-0.84) 


ASA 


12.26 (14.35) 


0.53 (-0.15) 


-1.31 (-1.48) 


-2.14 


(-0.98) 


Neutr. (GPM) 


5.49 


1.22 


0.01 


-0.73 






Connolly- Williams 


method 








ASA+M 


12.68 


1.31 


-0.02 


-0.73 




ASA+M (0-L) 


13.70 


0.49 


-0.86 


-1.39 




ASA 


14.33 


0.28 


-1.72 


-1.92 




US-PP 


12.81 


1.30 


0.69 


-0.40 






Direct calculations 


from (19) 








ASA+M 


12.45 


0.47 


-0.49 


-0.65 





TABLE V. Ordering energies of the equiatomic NiPt alloys obtained from the pair SGPM 
effective interactions. The values obtained in the direct calculations (see Table II) are given in 
parentheses. 



Structure 


ASA+M 


ASA 


Neutr. 


Ll 


-7.73 (-8.05) 


-7.46 (-8.92) 


-3.61 (-3.82) 


Z2 


9.22 ( 8.66) 


8.68 ( 9.85) 


4.19 ( 4.14) 


"CH" 


-7.15 (-6.86) 


-6.62 (-7.87) 


-1.77 (-1.67) 


Lli 


-2.96 (-2.38) 


-3.75 (-2.58) 


-3.09 (-3.54) 



TABLE VI. The on-site screening constant, ag Cr (0) in the first 3 layers of (111), (100) and 
(110) fee surfaces. 



facet 


A =1 


2 


3 


fcc(lll) 


0.805 


0.730 


0.728 


fcc(100) 


0.841 


0.732 


0.729 


fcc(110) 


0.840 


0.756 


0.732 
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TABLE VII. a4 A r '(i?i) at the first four coordination shells of the (111), (100) and (110) fee 



surfaces. 



A 


facet 


A' =1 




2 


3 


4 








Ri 








1 


fcc(lll) 


0.1160 




0.1262 


- 


- 




fcc(100) 


0.1007 




0.1241 


- 


- 




fcc(llO) 


0.1026 




0.1176 


0.1285 




2 


fcc(lll) 


0.1297 




0.1295 


0.1287 






fcc(100) 


0.1263 




0.1322 


0.1298 






fcc(llO) 






0.1325 


0.1306 


0.1318 


3 


£ ( 1 1 1 \ 

tcc(lll) 






0.1319 


0.1317 


0.1307 




fcc(100) 






0.1325 


0.1306 


0.1318 




fcc(llO) 


0.1290 




0.1297 


0.1285 


0.1312 








R 2 








1 


fcc(lll) 


- 




0.0012 


- 


- 




fcc(100) 


-0.0111 




- 


-0.0022 


- 




fcc(llO) 


-0.0111 






0.0040 




I 


/ 1 i 1 \ 

tcc(lll) 


-0.0035 






0.0012 






fcc(100) 






-0.0003 




0.0017 




fcc(llO) 






-0.0021 




-0.0011 


o 



£ ( 1 1 1 \ 

fcc(lll) 






-0.0003 




0.0001 




fcc(100) 


0.0027 






0.0016 






fcc(llO) 


-0.0039 






-0.0023 










Rs 








1 


fcc(lll) 


-0.0139 




-0.0099 


-0.01441 


- 




fcc(100) 


- 




-0.0080 


-0.0156 


- 




fcc(llO) 


-0.0131 




-0.0090 


-0.0101 


-0.0111 


2 


tcc(lll) 


n ai A/i 

-0.0104 




-0.0096 


n nan r 

-0.0095 


-0.0111 




fcc(100) 


-0.0097 






-0.0088 


-0.0096 




fcc(llO) 


-0.0079 




-0.0069 


-0.0107 


-0.0117 


3 


fcc(lll) 


-0.0099 




r\ AAA 'T 

-0.0097 


-0.0111 


A A1 1 O 

-0.0113 




fcc(lOO) 


-0.0091 




-0.0087 




-0.0103 




fcc(llO) 


-0.0124 




-0.0107 


-0.0118 


-0.0108 








R 4 








1 


fcc(lll) 


-0.0058 




- 


-0.0087 


- 




fcc(lOO) 


-0.0043 






-0.0103 






fcc(llO) 


-0.0031 






-0 0063 




2 


fcc(lll) 






-0.0031 




-0.0042 




fcc(lOO) 






-0.0033 




-0.0030 




fcc(llO) 






-0.0047 




-0.0036 


3 


fcc(lll) 


-0.0024 






-0.0034 






fcc(lOO) 


-0.0019 






-0.0030 






fcc(llO) 


-0.0029 






-0.0058 
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